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Abstract. We compute analytically the density e7V,M(A) of Schmidt eigenvalues, 
distributed according to a fixed-trace Wishart-Laguerre measure, and the average 
Renyi entropy (Sg) for reduced density matrices of entangled random pure states 
with orthogonal symmetry (/3 — 1). The results are valid for arbitrary dimensions 
N — 2k, M of the corresponding Hilbert space partitions, and are in excellent 
agreement with numerical simulations. 
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1. Introduction 

Entanglement is one of the most distinctive features of quantum systems. Recently it 
has attracted much attention in view of possible applications to quantum information 
and quantum computation problems [HE]- In these domains, one is often interested 
in creating states with large entanglement, thus raising the question how to give a 
quantitative measure of entanglement. Pure bipartite systems (defined below) constitute 
a typical example where well-behaved entanglement quantifiers can be defined, such as 
the von-Neumann or Renyi entropies of either subsystem [2], the so called concurrence 
for two-qubit systems [3] or other entanglement monotones [HE]. 

Introducing a source of randomness in quantum entanglement problems is the key 
to address typical properties of such states: in this paper we focus on random pure 
quantum states in bipartite systems, where many analytical results have been obtained 
in recent times (see e.g. [6| for an excellent review). 

More precisely, we consider a bipartite partition of a A^M-dimensional Hilbert space 
:K(^^) as :K(^*^) = :K^f ® :K^j^'\ where we assume without loss of generality that 
< M. For example, A may be taken as a given subsystem (say a set of spins) and 
B may represent the environment (e.g., a heat bath). Take a quantum state {ip) of 
the composite system and let {i"^) and \ct^) be two complete basis of JC^'' and "K^^^ 
respectively. The state {ip) can then be expanded as a linear combination 



1=1 a=l 

whose coefficients Xj^^'s form the entries of a rectangular (A^ x M) matrix X. 
Possible features of we are considering here are the following: 

• Entanglement: we say that {ip) is an entangled state if it cannot be expressed as 
a direct product of two states belonging to the two subsystems A and B. In other 
words, in order for \ip) to be fully unentangled, the coefficients Xi^a must have the 
product form Xi^a = ca^a for all i and a in a certain basis. In this case, the state 
\'tp) = 19"^) ^\9^) can be written as a direct product of two states 19"^) = ai\i^) 
and l^^'^) = Xlali ba\(^^) belonging respectively to the two subsystems A and B. 

• Randomness: suppose that the expansion coefficients Xt^a are random variables 
drawn from a certain probability distribution. In this case, we say that \ip) is a 
random state, and here we focus on the simplest and most common case where Xi^aS 
are independent and identically distributed (real or complex) Gaussian variables. 

• Purity: the density matrix of the composite system is simply given by p = 

with the constraint Tr[p] = 1, or equivalently {ipH') = 1. Note that the composite 
system may instead be in a statistically mixed state, with a density matrix of the 
form 



N 



M 




(1) 




(2) 



k 
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where \4'k)'s are the pure states of the composite system and < pk < I are the 
probabihties that the composite system is in the fc-th pure state, with YlkPk — 1- 
We will not consider this case here, and we refer to [7] and references therein for 
recent results on mixed states. 

Let \ip) be an entangled pure state of a bipartite quantum system. Its density 
matrix can then be straightforwardly expressed as 

where the Roman indices i and j run from 1 to and the Greek indices a and /3 run 
from 1 to M. We normalize the pure state to unity so that Tr[p] = 1. 

Tracing out the environmental degrees of freedom (i.e., those of subsystem B) leads 
to the definition of the reduced density matrix pA = Ttb[p]- 

M 

p^ = Trs[p] = 5^(«^|p|«^). (4) 

a=l 

Using the expansion in Eq. ([3]) one gets 

N M N 

= E = E W^^^''){A (5) 

where VFj/s are the entries of the N x N matrix W = XX^. In analogous way, one 
could obtain the reduced density matrix ps = Tryi[p] of the subsystem B in terms of 
the M X M matrix W = X+X and find that W and W share the same set of nonzero 
(positive) real eigenvalues {Ai, A2, . . . , Xn}- In the diagonal basis, one can express pA as 

N 

PA = l]A.|Af)(Af| (6) 

i=l 

where |Af)'s are the normalized eigenvectors of W = XX"'' and similarly for ps- The 
original composite state in this diagonal basis reads: 

N 
i=l 

Eq. is known as the Schmidt decomposition, and the normalization condition 
(iplip) = 1, or equivalently Tr[p] = 1, imposes a constraint on the eigenvalues. 

It is useful to remark that while each individual state \Xf-) ® |Af ) in the Schmidt 
decomposition in Eq. ([7]) is unentangled, their linear combination lip), in general, 
is entangled, and therefore the state cannot, in general, be written as a direct 
product \ip) = ® l^'^) of two states of the respective subsystems. Knowledge of 
the eigenvalues {Ai, A2, . . . , Xn} of the matrix W is essential in providing information 
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about how entangled a pure state is. Typical entanglement quantifiers include the Renyi 
entropy of order q > 1 



which converges to the von Neumann entropy Svn = —Yli=i'^i^^^i q ^ 1. The 
Renyi and von Neumann entropies attain their minimum value when one of the 
eigenvalues reaches its maximum value 1 and all the others are zero, which corresponds 
to completely unentangled states, while they attain their maximum value InA^ in the 
situation where all eigenvalues are equal (Aj = 1/N for all i). In this case, all the states 
in the Schmidt decomposition (JTj) are equally present and the state is maximally 
entangled. 

So far, we have considered an arbitrary pure state in Eq. ([T]) with fixed coefficient 
matrix X = [xj^Q,]. This state is called random if the coefficients are drawn from an 
underlying Gaussian distribution (real or complex) Prob[X] oc exp [-fTr(XtX)] where 
the Dyson index /3 = 1,2 corresponds respectively to real and complex X matrices. 
While generally {xi^a} are complex, real coefficients are important for systems enjoying 
a time-reversal (or any anti-unitary) symmetry. In these cases, it is known that the 
eigenfunctions can be chosen to be real, and the corresponding ensembles are the 
'orthogonal' ones (/3 = 1). Exact results for the statistics of random orthogonal states 
are very scarce [HI |9] . It is the goal of this paper to fill this gap and to present exact 
results for the average density of Schmidt eigenvalues (one-point function) and the 
average Renyi entropy, valid for arbitrary dimensions N = 2k, M of the corresponding 
Hilbert space partitions. 

Conversely, analytical results for spectral statistics of random pure states with 
broken time-reversal symmetry (/3 = 2) abound. The joint probability density (jpd) 
of Schmidt eigenvalues was derived by Lloyd and Pagels [lOj (see eq. f|T^ below), and 
using this result Page [H] computed the average von Neumann entropy for large N, M 
and found: 



Since In(A^) is the maximal possible value of von Neumann entropy for the subsystem 
A, in the limit when M ^ N, the average entanglement entropy of a random pure state 
is close to maximaiji. Later, the same result was shown to hold for the /3 = 1 case |12j . 
In the same paper. Page also conjectured from numerical experiments that the average 
von Neumann entropy for finite N, M and /3 = 2 should read 




(8) 



(§vn) ~ In(iV) 



N 



(9) 



2M 




(10) 



k=N+l 



I Note, however, that the probabihty of the maximally entangled microscopic state (where all Schmidt 
eigenvalues are close to each other) decays very quickly as N increases, a result that is based on the 
exact evaluation of the full large deviation tails [22] . 
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a result that was independently proven by many researchers soon after [13] also in 
a non-extensive setting [H]. Recently, many efforts have been directed towards the 
study of other statistical quantities for finite {N, M), and full distributions of interesting 
observables. We mention here: 

• Density of Schmidt eigenvalues (one-point function) for (3 = 2 and finite {N,M), 
derived independently in [T5] and [T6] : 

• Universality of eigenvalue correlations for /3 = 2 [T7j; 

• Distribution of minimum eigenvalue for /3 = 1,2 and finite {N,M), derived in [8] 
where a conjecture by Znidaric fTSj was proven (see also |9j for a related result); 

• Average fidelity between quantum states [19] and distribution of so-called G- 
concurrence ^ for /3 = 2; 

• Distribution of so-called purity (i.e. §2) for small [20], and phase transitions in 
its Laplace transform for large N [21] ; 

• Full distribution of Renyi entropies (including large deviation tails), computed in 
[22] for large N = M and all (3s using a Coulomb gas method. As a byproduct, 
the authors also obtain in [22] the average and variance of Renyi entropy valid for 
large A^ = M as @: 

(S,) ^InAT-i^n^lnAT-i (11) 
Var(S,) « ^ (12) 

where: 

4''r(^ + 1/2) 
= v^r(, + 2) 

We will compare in Section |3] the asymptotic result ffTTl) with our exact formula for 
the average (S^) for /3 = 1 (see eq. P2l) ) and find that fH2l) converges to (|TT1) very 
quickly for low q, thus including the most relevant cases q = 1,2. Conversely, the rate 
of convergence progressively deteriorates as q increases (see Section [3]). 

In order to proceed, we now summarize the basic ingredients of the calculation. 
The joint distribution of Schmidt eigenvalues Aj G [0, 1] [10] reads: 

. . . , A.) = CZ,S E^' - 1 n A? n - A.I* (14) 

\j=l / i=l j<k 

where C^]y^ is a normalization constant known exactly for any /3, and the Dyson index 
P = 1,2 identifies respectively systems with preserved (orthogonal) or broken (unitary) 
time- reversal symmetry. The delta function guarantees that Tr[p] = 1 and implies that 
the typical eigenvalue scales as A ~ 1/A^. 

Another jpd of eigenvalues which is closely related to ( |T4l) is from the Wishart- 
Laguerre ensemble of random matrices [23l [23] of the form W = X^X, where X is a 

§ Note that the limit g — >■ 1 of eq. (fTTjl is consistent with eq. ([9]) already derived by Page. 
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Gaussian M x N matrix with real or complex entries. The joint distribution of the 
nonnegative eigenvalues of W is known [25J 

a'l"')(A„...,A„) = |0-e-t^".^- nAp"-"'-' ni^i-A*!" (15) 

1=1 j<k 

where k]^\j is a known normalization constant. Therefore, the jpd (fT4l) can be seen 
as a fixed-trace version of the Wishart-Laguerre ensemble. The presence of a fixed- 
trace constraint has crucial consequences on the spectral properties of random matrix 
ensembles [2^ [TT] . The goal of this paper is to compute exactly the one-point marginal 
(average density) Qnm^^) for /3 = 1, defined as: 

1 ^ 

Qn,m{x) = (^—^5{x-\)^ (16) 

i=l 

where the average (■) is taken with respect to the measure (fT4|) . Writing down this 
average explicitly, one is led to: 

^?iv,M(Al) := di^u / d\2--- d\^b V A, - 1 n Af^ n |A, - Afcl (17) 

where v = M — N . Computing this (A^ — l)-fold integral is the main technical challenge. 
Note that in the large N,M limit with c = N/M fixed, the average density can be 
computed for all /3s using a Coulomb gas technique [22] and has the scahng form: 

Qcix) = Ng*^{Nx) (18) 

where: 

Qtix) ■■= ^^(L(+)(c)-x)(x-L(_)(c)) (19) 

where the edge points L(^±) = ci^^/xjc ± 1)^. 

The average density is important in order to obtain averages of so-called linear 
statistic^ with a simple one- dimensional integration as: 

N „i 

Y,f{\))=N \ d\ q^,m{\)!{\) (20) 

In particular, we consider f[x) = x'^ for the Renyi entropy in Section [31 

The paper is organized as follows. In Section [21 we compute the one-point density 
for /3 = 1 and any = 2k, M using a Laplace transform method, which is technically 
transparent and avoids the unnecessarily heavy formalism used for earlier derivations 
of the (3 = 2 case [16] . In section [31 we use the obtained result to compute the average 
Renyi entropies for /3 = 1 and compare them with the asymptotic formula ( |TT|) for large 
N . Then we provide some conclusions in Section [H 

II A linear statistics is a quantity of the form = 'Yl!i=i fi^i): where f{x) is any smooth function. 
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2. Density of Schmidt eigenvalues for the orthogonal case (3 = 1 (N even) 

In this case the jpd ( HM of Schmidt eigenvalues reads: 

m^J^ni^^-^^i (21) 

i=l / i=l j<k 

where we put Cn,m = C^nm ^"^^ v = M — N . 

The goal is to compute the density of eigenvalues for finite (A^, M), i.e. the marginal 



Qn,m (Ai) := 



/ dA2 ■ ■ ■ d\j,6 V A, - 1 n Af^ n |A, - Afcl (22) 

which is normalized to 1, i.e. Jq dzg^^Mi^) = 1- In the orthogonal case the 
normalization constant reads Cn m = ,-,^-1 ^(^^/^)(v^/^) 



n iLV r((M-j)/2)r(i+(iv-j)/2) • 
We first define Qn,m{.^i) = ^Ar,Af(Ai, 1), where 



r ( ^ A ^ .-1 

^iv,M(Ai, t) := Civ,A/ / d\<,--- dXNS V A, - t n A7^ TT |A, - Afel (23) 



is an auxiliary function that we are going to compute exactly. 
We next take the Laplace transform of (!23|) : 

poo n ^ ^„j^ 

/ dtQN,M{Kty-'' = Cn,m dX2...d\Ne-'^^^^^Wx;^W\\-\u\ (24) 

where in the r.h.s. we have extended the range of integration to the full positive semiaxis. 
This is harmless in view of the unit norm constraint. The integral on the r.h.s. can be 
written in the form: 

/ . . . dA„.-££.^. n xp" n - = ^f,,S:M./2 A's'(2^^o (25) 

where: 



= 2-/^(^/^/2)- nr (^) r (^l^jy±A±l) (26) 

is the normalization constant of the jpd of eigenvalues of a Wishart-Laguerre (WL) 
ensemble with /3 = 1 (which can be derived from the Laguerre-Selberg integral (see e.g. 
[28]), and 

r ^ -1 

q7m{>^i) = N{Kn,m)-' / dX2...dXNe-"^^^^^^T\Xr^T\\\,-Xu\ (27) 

>/[o,-)^-^ t=l J<k 

is the one-point density of the WL ensemble, normalized to A^. 
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The spectral density (one-point function) for the WL ensemble Q^^MiXi) at even 
is known: 

= Ma;(-i)/2e-/2 r dx'sgn{x - x'){x'Y''-'^/^e~^'/^S{x,x';iy,N) (28) 
4 Jq 

where: 



S{x, x'- V, N) := J2 77T^{^.>i(^')^.'(^) - ^ • (^')} (29) 



i=o ^-^ 

where L'^{x) are Laguerre functions defined by the sum: 



iV! Z^A:!r(z/ + /t + l)'^ ^^"^ 

(where (x)„ = r(x + n)/T{x)) and sgn(z) = z/\z\. The explicit formula (128|) can be 
most conveniently derived by taking the yU — > limit of eq. 4.14 in [2S]. Equivalent 
but less handy expressions can be found in [30], while the general formalism based on 
skew-orthogonal polynomials is in |2H |3T] . 

In order to take the inverse Laplace transform of ( l25l) . some work is needed. First, 
we make a change of variable x' = xz in f l28|) . obtaining: 

g'-;^l^{x) = sgn{l - z)z^''-''>/^e-^'/^S{x,xz;u,N) (31) 



4 

For later convenience, we now define and compute the following inverse Laplace 
transform: 

^k{t,x,z;N,,N2;u) := L'^ [sV^^(i+^)L^^(2sx)L^^(2sx^)] (t) (32) 

Using the general definition of Laguerre functions fl5U]) and the following elementary 
Laplace inverse: 

.-.|„«e-„*, ^ "-"'^'liy-^' (33) 

(where 6{x) is the Heaviside step function), it is straightforward to get: 

, r(i^ + A^i + i)r(z/ + A/'2 + 1) 

^y^y^ {-N,U-N2h ^ 

^^j^mU\T{u + m + l)T{u + i + l)T{-k - m - £) 

{2xT^'z' it - x(l + 2))-'-'"^'-' ^ (t - a;(l + z)) (34) 

Combining everything together, we obtain for ^Ar,M(Ai,t): 



0N,M (Ai,t) 



2l+AfAr/2JY 



X [^,(t, Ai, 2;; J, J + 1; I/) - ^.(t, Ai, z;j + 1, j; z/)] (35) 

where: 

K = M -N- MN/2 + 1 (36) 
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The sought density {t = 1) can then be written in the compact form: 

N-2 j j+1 

qnM^) = ^N,M E E E ^^ey^'''^' X 

j=0 m=0 e=o 



X dz sgn(l - z) z^"-^^^'' (1 - xil + z))-"-'"-^-' (l)iUz)Oil - + z)) 



where: 

_ Cn,mKn,M , . 

u) ^ r(^ + j + 2)2-n-jU(-j-i)^ .3g^ 

^"^ j\mW.T{u + m+l)T{iy + e + l)T{-K-m-i) ^ ' 

ct>,^{z) = /-z^ (39) 

The integral in z can now be performed exactly. Let 

E{a,b]x)= dzsgn{l-z)z''{l-x{l + z))'', Re[a, 6] > -1(40) 

Jo 



One has: 

{B(a + 1,6 + 1), forl/2<x<l 
-B(a + 1, 6 + 1) + 2B(x/(l - x), a + 1, 6 + 1), 
for < X < 1/2 

where B(a,6) = r(a)r(6)/r(a + b) is Euler's Beta function and B{z,a,b) = 
du u°'^^{l — n)^~^ is the incomplete Beta function. 
Eventually one gets for the density of Schmidt eigenvalues for /3 = 1: 

N-2 j i+1 r / _ 1 \ 

f?7V,M(x) = J^NM cJ-^'"""^' - ( ^ m - £ - 1; X j 



j=0 m=0 e=o 



+ m, — K — m — £ — 1: X 



< X < 1 (41) 



2 

Equation (14T!) is the main result of this sectioii^. The obtained exact formula is the 
starting point to compute averages of linear statistics using formula (!20l) . In the next 
section, we are going to compute the average Renyi entropy at finite A^, M and compare 
it with the exact asymptotic result for large N = M obtained in [22]. Note that 
(A) = Jq dX XqnmW — agreement with the general scaling argument that 

typically A ~ 1/A^ due to the trace constraint J2i=i -^^ — ^■ 

In fig. [Dwe plot the density dS]) for /3 = 1, 2 for = M = 6, and in fig. |2]the 
density f HT]) for = 6,M = 12 together with the large N density f|T8|) . In fig. [3] we 
compare the theoretical density with numerical results, obtained as follows [32| [33]: 



% Following analogous but much quicker steps, one can also derive the already known one-point density 
for /3 — 2 in a much simpler way. For example, for N — M we obtain 

(/^=2), , T{N^) ^ i-k), i^kU ^^+"(1 - x)^''-^-^-"^ 
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(i) we generate n = 5 ■ 10^ real Gaussian M x N matrices X (where = 6, M = 8). 

(ii) for each instance we construct the Wishart matrix W = X^X. 

(iii) we diagonahze W and collect its real and non-negative eigenvalues {Ai, . . . , Xn}- 

(iv) we define a new set of variables < Aj < 1 as Aj = Aj/ J2iLi Aj, for i = 1, . . . , iV. 
The set of variables Aj is guaranteed to be sampled according to the measure f|T^ . 

(v) we construct a normalized histogram of Aj. 

The agreement between theory and simulations is excellent. 

Pn,mOO 
12 

10 

8 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 1. Density of Sclimidt eigenvalues for 13 — 1 (solid blue line, eq. (|41l) ) and 
/3 = 2 (dashed red line), both for N — M — 6. 



3. Average Renyi entropy 

The average Renyi entropy (S^) is computed from the density f HT]) as: 



1 



l-q 



log 



N I dx x'^qn m{x) 
' ^ ' 



The integral d{q) can be computed easily. First, define the function: 

(p{a,b,c,d) := f dz z^'il + zfB{z,c,d) 
Jo 

Then the sought formula for J(g) reads: 

N-2 j j+1 

d{q) = y^NM E E E [^^(^ + ^) - ^™("^ + ^)] 

j=0 m=0 e=o 



(42) 



(43) 



(44) 
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Pn,mM 

6r 




0.6 



0.8 



1.0 



Figure 2. Density of Sclimidt eigenvalues for /3 = 1 (eq. (jUJ) for N — 6,M ~ 12 
(solid line). In dashed red, the corresponding large N density ([TS]) for c — N/M = 1/2. 




Figure 3. Density of Schmidt eigenvalues for /3 = 1 (eq. (|4T|) . solid line) for 
N = Q, M — 8 compared with numerical simulations (red dots). 



where: 



eria,q) :=B 



u - 1 



r + 1, —K — a ) Gr{a, q) + 



V + 1 V -I Z/-1 
+ 2^9 g H h « — 1 — r, — h r — k — a, — h r + 1, — k — a 



Gr(a;, q) ■.= B[q + a-r + + r- K- a + l 
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-U+i 2F1 iq + a - r -\ q + iy-K + l;q + a-r + 



2 



Here, 2-^1 (a, b; c; x) is a hypergeometric function defined by the series: 



00 



2Fi(a, 6; c; x) = ^ 



(45) 



In fig. m and [5] we compare respectively (§2) (average purity) and (§6o) as a function of 
A^(= M) for /3 = 1, 2 with the large N asymptotic formula (fTTj) from [22]. 
We find that: 

(i) The average Renyi entropy for systems with time- reversal symmetry (/3 = 1) is 
always lower than systems of the same size where this symmetry is broken (/3 = 2). 
This fact is in agreement with recent findings [31] about so-called 'single-particle' or 
one-magnon states, where real states have lower entanglement measured in terms 
of two-spin entanglement content than the case of complex states. We have checked 
that this feature persists for N 7^ M, where a large N formula for general q is not 
yet available. 

(ii) For low q, the finite and large results are in excellent agreement already for 

~ 6, for both /3 = 1,2. This means that for the most relevant cases of average 
von Neumann entropy (g — t- 1) and purity (g = 2), one can safely use eq. (fTTIl 
from [22] as an excellent approximation for any practical purposes. Conversely, the 
quality of the approximation decays as q increases, up to the limit Sqo — — In Amax 
(where Amax is the largest Schmidt eigenvalue), and one has to consider larger and 
larger subsystems in order to reach a satisfactory agreement (see fig. [5]). The 
discrepancy between the /3 = 1 and /3 = 2 is also more pronounced in the case of 
high q, and the convergence to the asymptotic limit is much slower for the (3 = 1 
case. 

4. Conclusions 

In summary, we have computed exactly the density ^Ar,Af(A) of Schmidt eigenvalues for 
bipartite entanglement of random pure states with orthogonal (time-reversal) symmetry 
(/3 = 1). The result is valid for any finite dimensions N < M (with N even) of the 
corresponding Hilbert space partitions. Using the exact formula we derived and a simple 
linear statistics, we compute the average Renyi entropy (Sg) for the /3 = 1 case, which 
was previously unavailable. We find that the exact values for the averages at = M 
converge very quickly to the asymptotic N ^ 00 formula derived in [22j for low values 
of the parameter q, thus including the most relevant cases of the von Neumann entropy 
[q — 7- 1) and the so-called purity {q = 2). As q is increased, the speed of convergence 
deteriorates for both /3 = 2 and /3 = 1, and the latter value for the average Renyi entropy 
is consistently lower than the former for the same values of parameters q, N, M, even at 



N ^M. 
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1.0 - 



0.5 - 




Figure 4. Average purity (§2) as a function of N{= M) for /3 = 1 (green dots) and 
(3 = 2 (red dots), compared with the exact asymptotic result eq. pT|) for — > 00 [22] , 



<S60> 




Figure 5. Average Renyi entropy (§6o) (<Z ~ 60) as a function of A^(= M) for /3 = 1 
(green dots) and (3 = 2 (red dots), compared with the exact asymptotic result eq. (fTT|) 
for N ^00 gl]. 
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